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Abstract. 

We discuss the present bounds on the relativistic energy density in the Universe 
parameterized in terms of the effective number of neutrinos N^ft using the most recent 
cosmological data on Cosmic Microwave Background (CMB) temperature anisotropies and 
polarization, Large Scale galaxy clustering from the Sloan Digital Sky Survey (SDSS) and 2dF, 
luminosity distances of type la Supernovae, Lyman-a absorption clouds (Ly-a), the Baryonic 
Acoustic Oscillations (BAO) detected in the Luminous Red Galaxies of the SDSS and finally, 
Big Bang Nucleosynthesis (BBN) predictions for 4 He and Deuterium abundances. We find 
N£" = 5.2^2 2 fr° m CMB and Large Scale Structure data, while adding Ly-a and BAO 
we obtain N^P = 4. 6+^5 at 95 % c.L. These results show some tension with the standard 
value Nf = 3.046 as well as with the BBN range Nf = 3.1±^ at 95 % c.L, though the 
discrepancy is slightly below the 2-er level. In general, considering a smaller set of data weakens 
the constraints on N£" . We emphasize the impact of an improved upper limit (or measurement) 
of the primordial value of 3 He abundance in clarifying the issue of whether the value of Nf^ 
at early (BBN) and more recent epochs coincide. 
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1. Introduction 

The energy density in relativistic species pr is one of the key parameters defining a specific 
cosmological model and is usually parameterized in terms of the effective neutrino number 



with p 7 (T) the photon energy density and T v the neutrino temperature. There is a very 
neat result for its value after e + — e~ annihilation, namely N** = 3.046 assuming standard 
electroweak interactions, three active neutrinos and including the (small) effect of neutrino 
flavor oscillations [1] . On the other hand, Nj™ which in complete generality is a function of the 
particular stage of the expanding Universe, i.e. of the photon temperature, can be indirectly 
constrained by several cosmological observables. The final yields of both 4 He and D produced 
during Big Bang Nucleosynthesis (BBN) are sensitive probes of the total amount of radiation 
energy density which drives the expansion rate via Friedmann equation, as well as of the 
neutrino phase space distribution which enters the weak rates governing the eventual neutron 
to proton number density ratio. At later stages, the imprint of different values of N£" can be 
extracted from the temperature and polarization maps of photons during the formation at the 
last scattering surface of the Cosmic Microwave Background (CMB) or the power spectrum of 
Large Scale Structures (LSS), mainly because its role in shifting the matter-radiation equality 
and changing the early Integrated Sachs- Wolfe effect. For recent reviews on these issues see 
e.g. [2J and [3]. 

The huge amount of observations from CMB satellite, balloon-borne and ground based 
experiments, pE]-[5], galaxy redshift surveys, PUD], distant type la Supernovae [IT] , to only 
mention those of larger impact, is pointing out a rather consistent standard Cosmological Model 
(see however, [12] for a different viewpoint). On the other hand, all these data can be fruitfully 
used to constrain non standard physics at the fundamental level, such as classes of elementary 
particle models predicting a different radiation content in the Universe. 

Indeed, there are several theoretical models which have been considered to account for 
different values of Nfp at different stages as the physics involved at the BBN and at the CMB 
epochs could be different [131 Eh see e -g- [E] f° r the baryon energy density values inferred 
from BBN and CMB data sets separately. Many of these models involve additional sterile 
neutrino species which either are generated in presence of a large lepton number and their 
production is therefore suppressed avoiding constraints from CMB and LSS, or are produced in 
the decay of an active neutrino together with a majoron [16] . These models, in general, assume 
a (light) massive sterile neutrino, motivated by the LSND oscillation claim [17J. However, 
current cosmological data excludes at a high significant level the sterile neutrino hypothesis as 
an explanation for the LSND signal [18]. Nevertheless, it remains still viable to enlarge the 
neutrino sector with massless sterile species by imposing a discrete symmetry which prevents 
the extra sterile neutrinos to acquire masses. For example, it has been recently suggested in [19] 
a Leptocratic model, which enlarges the Standard Model by including a non-anomalous U(l) v 
symmetry and predicting the existence of light (m„ <C 1 eV) quasi-sterile neutrinos. Indeed, 
from a theoretical point of view, all the models which contain massless or very light sterile 
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neutrinos pose a major challenge, namely to understand the sterile neutrino lightness, despite 
being a Standard Model singlet. An interesting possibility is the model discussed in [20], with 
neutrinos acquiring a mass from a broken lepton flavor symmetry via operators of the form 
O = £Nh(cj)/M) p , where £ is the left-handed lepton doublet, h is the Standard Model Higgs 
field, M a mass scale larger that the electroweak scale, p a positive integer and finally, N is a 
right handed neutrino. Such a model contains a number of pseudo-Goldstone bosons (0) which 
could modify the relativistic energy density at the eV era (CMB) while keeping BBN essentially 
unchanged if these extra particles only interact with Standard Model particles via operators of 
the form O. This model have predictions for Nf® at the CMB epoch in the range 3 4- 4. For 
a recent combined study constraining the Yukawa couplings of neutrino to low mass (pseudo) 
scalars, see [21]. 

Another possible scenario is to allow for a violation of the spin-statistics theorem in the 
neutrino sector, so that the neutrino distribution function interpolates between the fermionic 
and the bosonic results, increasing the number of degrees of freedom for CMB from up to 
3/7, due to the presence of bosonic neutrinos [22J. 

A different possibility is finally, to consider either a quintessence field with a tracking 
potential behavior or an extra interaction among the dark energy and radiation (or dark matter) 
sectors or even a Brans-Dicke field with non trivial potential [23] which could mimic the effect 
of adding extra radiation between the BBN and CMB epochs. Indeed, we presently know very 
little about the role of any possible dark energy between BBN and CMB eras, despite the fact 
that the amount of dark energy at both epochs is very well constrained [24|. 

These numerous theoretical studies, also motivated by the improvement of the data 
accuracy which is expected from the next generation of planned experiments like PLANCK, shows 
the fact that there is quite an intense research activity which intrigues several scholars and try 
to answer the following two main questions 

i) what is the value of which is favored by cosmological observations with respect to the 
standard result? 

ii) is there any hint for a time dependent value of Nf™ , which takes different values at, say, 
the BBN epoch (photon temperature of order MeV) and CMB last scattering (1 eV)? 

In this paper we address both these issues. Though several analysis of the kind are already 
present in the literature [13], [25]- [38], we provide updated results using a broad set of the most 
recent available data from CMB, LSS, as well as using Baryon Acoustic Oscillations (BAO) 
detection of [39], small scale spectrum from Lyman-a forest [401 |4T] and finally, reconsidering 
the constraints from BBN in view of a new Deuterium measurement reported in [42] and using 
an improved BBN numerical code described in [2] . In particular, our study has been stimulated 
by a a recent paper [43] where a mild disagreement (at the 2.4-cr level) is found between the 
(large) value of at late epochs (CMB and LSS) and both the standard value and that 
singled out by BBN. Section II is devoted to a summary of the data analysis method, while in 
Section III we discuss our results. Finally, we draw our conclusions and outlooks in Section IV. 
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2. Data analysis method 

The method we adopt is based on the publicly available Markov Chain Monte Carlo package 
cosmomc jH] with a convergence diagnostics done through the Gelman and Rubin statistic. 
We sample the following seven- dimensional set of cosmological parameters, adopting flat priors 
on them: the baryon and Cold Dark Matter densities, Ub = VL^h 2 and uj c = Q c h 2 , the ratio 
of the sound horizon to the angular diameter distance at decoupling, 9 S , the scalar spectral 
index ns, the overall normalization of the spectrum A at k = 0.05 Mpc -1 , the optical depth to 
reionization, r and the effective number of massless neutrinos N™. 

We consider two sets of cosmological data. In the first set, we include the three-year 
WMAP data [4] (temperature and polarization) with the routine for computing the likelihood 
supplied by the WMAP team. Together with the WMAP data we also consider the small- 
scale CMB measurements of CBI [5], VSA [6J, ACBAR [7J and BOOMERANG-2k2 [8J. In 
addition to the CMB data, we include the constraints on the real-space power spectrum of 
galaxies from the SLOAN galaxy redshift survey (SDSS) [9J and 2dF [TU] . and the Supernovae 
Legacy Survey data from [TTj . Finally we include a prior on the Hubble parameter from the 
Hubble Space Telescope Key project [45] ■ In the second dataset we add to the data quoted 
above the constraints from the Baryonic Acoustic Oscillations (BAO) detected in the Luminous 
Red Galaxies (LRG) sample of the SDSS [39] and we include measurements of the small scale 
primordial spectrum from Lyman-a forest clouds [4U| W\\. 

The first dataset is based on well-established data while the second include new results 
that, especially in the case of Lyman-a data, should still be confirmed by further experimental 
data. Therefore, even if more stringent the results obtained from the second dataset should 
be considered as less conservative. In the framework of the cosmological model and using the 
data described above we produce constraints on Uf, and Njf? , after marginalization over the 
remaining 'nuisance' parameters. 

Constraints on these two parameters are also computed using standard BBN theoretical 
prediction as provided by the new numerical code described in [2], which includes a detailed 
treatment of the three active neutrino decoupling stage, as well as a full updating of all rates 
entering the nuclear chain based on the most recent experimental results on nuclear cross 
sections. We adopt for the D/H abundance ratio the new average of [42] obtained including 
a new measurement in a metal poor damped Lyman-a system along the line of sight of QSO 
SDSS1558-0031 

D/H = (2.82l°J 5 7 ) • 10- 5 , (2) 

We use the uncertainty as quoted in [42] , computed by a jackknife analysis. Indeed, the six 
individual D/H measurements used in their study show a scatter around the weighted mean. 
This might be a real effect due to non homogeneous D/H abundance, or rather the effect of 
underestimated errors in measurements. 

Finally, we also use the conservative range for the 4 He mass fraction Y p obtained in [46] 

Y p = 0.249 ± 0.009 . (3) 

The likelihood function is then constructed as in [2], by also taking into account the (small) 
D - 4 He correlation term in the inverse covariance matrix due to the effect of the propagated 
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Figure 1. Analysis from the first (more conservative) set of cosmological data (Left Panel) 
and including Lyman-a and BAO data (Right Panel). We show the marginalized contours at 
68% and 95% c.l. on the uib-Nf plane along with the analogous contours from BBN using D 
and 4 He (cnservative) experimental results reported in Eqs. @ and (dotted lines). 



uncertainty of the several nuclear rates on the theoretical nuclide abundances. 
3. Results 

In Figure [I] we show the marginalized constraints at 68% and 95% c.l. in the u^-Nf plane 
obtained using the first dataset (left panel) and the results from the second datasets which 
includes Lyman-a data and BAO (right panel). In both cases we also show the analogous 
contours from BBN using the experimental determination of 4 He and D discussed in the 
previous Section. As we can clearly see, an extra-relativistic component (Nf > 3) is favored 
by both sets of cosmological data. More specifically, we obtain the following marginalized 
constraints: u b = 0.0224 ± 0.0012 and Nf = h.2t 2 2 l at 95% c.l. from the first dataset and 
u b = 0.0224±0.0011 and Nf = 4.6±i;| at 95% c.l. when Lyman-a and BAO data are included. 
These data seems to indicate the presence of an extra-background of relativistic particles at 
2 — a level, in agreement with the previous analysis of |43j even if with a slightly lower statistical 
significance. The small difference from these previous results is explained by slightly different 
treatment of the data. We also note that [43J did not use the Hubble Space Telescope prior by 
default. 

In general, considering a smaller set of data would weaken the constraints on Nf . For 
example, combining only the WMAP data with the recent data from the analysis of LRG in the 
SDSS [U] gives Nf = 3.7±i$. Moreover, using different smaller datasets may give different 
bounds on Nf . This is shown in Tabled! Notice that combining WMAP data with the SDSS 
power spectrum prefers a large value for Nf ~ 8 though with a large uncertainty, while one 
gets a scenario compatible with the standard result Nf ~ 3 if the 2dF data are considered, as 
already noticed in previous analyses [43] . 

On the other hand, BBN is in very good agreement with the standard result Nf ~ 3. 
Using the D and 4 He mass fraction of Eqs. ([2]) and ([3]), after marginalization we find 
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Table 1. Constraints on N^p from different datasets. For each dataset we also show the 
corresponding best fit A% 2 between model with N^" — 3 and allowing for variation in Nfp . 
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u b = 0.0215ig;go29 and N f = 3.1±i:i at 95% c.l.. 

Very recently, a new analysis of 4 He mass fraction has been performed using a large and 
homogeneous sample of extragalactic H II region 08] and taking into account all main source 
of systematics via a Monte - Carlo approach. This results in a much smaller uncertainty than 
in 06]. Indeed, using two different He I line emissivities, they find Y p = 0.2472 ± 0.0012 
and Y p = 0.2516 ± 0.0011, respectively. These values when combined with the result on 
Deuterium of Eq. ([2]) lead to a much more constrained bound on the effective neutrino 
number, = 3.0 ± 0.2 and = 3.3 ± 0.2 at 95 c.l.. Marginalization over gives 

a basically unchanged result for the baryon fraction, Uf, = O^lltg^o an d = 0.0217io;oo2i- 
All constraints from BBN are also summarized in Table [U 

Notice that all cosmological data, spanning several decades in redshift or photon 
temperature provides results which are in agreement at the 2 — a level, with the only exception 
of the WMAP±SDSS dataset which singles out a very large value of , see Table [TJ We 
conclude that presently there is not a strong evidence for different values of N£" at the BBN 
and late epochs and that a slightly non standard value N£" ~ 4 seems to be preferred when 
combining all cosmological observables (first two rows of Table [1]) . 

It is also interesting to present our results in a slightly different way, by considering for 
each nuclide abundance its dependence on Nff . Assuming that the value of does not 
change between the BBN and the late CMB and LSS formation epochs we can compute using 
our numerical BBN code what are the values of Y p and D which are preferred at 68 % and 
95 % c.l. from the results obtained for Uf, and N£" from the two datasets of the previous 
Section. This is shown in Figure EJ along with the corresponding experimental values. Notice 
that as well known, both nuclei show a monotonically increasing behavior as function of the 
relativistic energy density. As from Figure dj we see that the measured values of Y p and D /H 
are in agreement with the cosmological constraints and a better agreement is obtained for the 
non standard value ~ 4. In particular we can derive the following marginalized ranges 
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Figure 2. The 4 He mass fraction (Top left Panel) and D/H ratio (Bottom left Panel) from 
the first (more conservative) set of cosmological data versus N^p (see text). We show the 
marginalized constraints at 68% and 95% c.L, while the regions between the two dashed 
horizontal lines are the 1-a experimental measurement bands. Plots on the right show the 
same results but adding BAO and Lyman-a data 



Y p = 0.273l{j:o25 and D/H= (3.4lS;jj) • 1(T 5 at 95% c.L for the first dataset and Y p = 0.2661^020 
and D/H= (3.1±^;f ) • lO" 5 in case we include in the analysis both BAO and Lyman-a data. As we 
mentioned, these results are only meaningful under the assumption that Nff does not change 
from the BBN till very recent epochs. A different conclusion would be obtained if we reason in 
the framework of models where Nfp can change in time as those summarized in Section I. In 
this case, data shows that nuclei abundances are best fit by a standard value N£" ~ 3 at the 
BBN epoch, while CMB and LSS prefer larger values due to some extra relativistic degrees of 
freedom, though with quite a large error. 

Besides 4 He and D, the most abundant nuclides produced during BBN are 3 He and 7 Li. 
As for the latter, there are two recent measurements, namely 

7 Li/H= (1.23 + 0.06) • 1(T 10 , (4) 

from [49] and 

7 Li/H = (1.26^) • 10' 10 . (5) 

from [50] • Both these values are a factor larger than three smaller than what is expected from 
a standard BBN scenario with only three active neutrinos contributing to NjjP , which gives 
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Figure 3. The 7 Li/H (Top Left Panel) and 3 He/H (Bottom Left Panel) from the first (more 
conservative) set of cosmological data versus . We show the marginalized constraints at 
68% and 95% c.L. We also show the 1-a experimental upper limit on 3 Hc/H of Eq. ©. The 
experimental results on 7 Li/H are out of the chosen range. Plots on the right show the same 
results but adding BAO and Lyman-a data 



7 Li/H=(4.7 ± 0.4) • 10~ 10 for u>b = 0.0224, where the uncertainty is due to propagation of 
uncertainties on all nuclear rates [2]. It is quite difficult that the theoretical estimate may be 
reconciled with the experimental value, even allowing for unlikely large values of N£" which 
would be in striking disagreement with other nuclei abundances. This is shown in Figure [3] 
(Top Panels) where we report the 7 Li abundance versus for the two datasets considered 
here. Even if N** is of order of 8 - 9 which is completely excluded by D and 4 He measurements, 
the predicted value for this nuclide is too high, suggesting a strong depletion mechanism of the 
Spite plateau value [50]. This is also suggested by the fact that a positive dependence on the 
metallicity has been found in [49], while there is only a weak evidence for such a dependence 
in the star sample studied in [50]. In any case, it is quite problematic to compare BBN 
predictions with 7 Li data, and further studies seem necessary in understanding possible source 
of systematics in POP II stellar modelling. 

The results for 3 He using the two datasets are shown in the bottom panels of Figure [31 
We recall that presently only an upper bound has been found for the primordial abundance of 
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this nuclide |5T] [jj] 

3 He/H < (1.1 ±0.2) • 1(T 5 , 



(6) 



since the experimental study is highly complicated by the fact that the primordial 3 He is 
strongly changed by stellar activity |52j. It is interesting to note however, that this upper 
bound is very close to the expected primordial value for a standard BBN scenario. Using 
u>b = 0.0224 and three standard neutrinos contributing to the relativistic energy density we 
get in fact from the BBN numerical code used in our analysis 3 He/H=(1.02 ± 0.03) • 10~ 5 , 
where again the uncertainty is due to the errors in the (experimental) values of the nuclear 
rates, largely dominated by the uncertainty in the 3 He destruction process 3 He(d,p) 4 He [2J. By 
looking at Figure [3] we see that any further improvement in the upper limit or an experimental 
estimate for its primordial abundance would have a big impact in constraining the value of Nf™ 
at the BBN epoch. For example, if the upper limit would decrease by 20 % this may represent 
a robust evidence in favor of a standard value N£" ~ 3 at BBN, so that any indication of larger 
values for this parameter by other cosmological observables would be a strong indication of a 
late production of relativistic degrees of freedom after BBN. From our analysis we also find 
that combining all data used so far we can derive the following conservative lower bound at 95 



where we have also included the effect of the nuclear rate uncertainty in the theoretical estimate. 
4. Conclusions and Outlooks 

In this paper we have presented an updated analysis of the constraints on the relativistic energy 
density in the Universe combining several cosmological data which measure the value of the 
effective number of neutrinos Nffi at different epochs in the evolution of the Universe. Nuclide 
abundances produced during BBN can tell us in fact how much energy was in the form of 
relativistic species when the photon temperature was in the range 0.01-1 MeV, while CMB and 
LSS data provide information at much later epochs. In the standard scenario after e + — e~ 
annihilation only three chiral active neutrinos contribute to , and any observed deviation 
from the value N** = 3.046 predicted in the framework of the Standard Model interactions 
would be a signal of new physics. 

Using WMAP and small scale CMB measurements, LSS and type la Supernova Surveys we 
find that N£" ~ 3 is allowed at 2-a though the data prefer somehow larger values. In particular 
we obtain Nfp = 5.2l|; 2 at 95 % c.L. Including Lyman-a and BAO data do not substantially 
change the result, Nff = 4.6l{;5 a ^ 95 % c.L. On the other hand, data on primordial D and 
4 He are in excellent agreement with ~ 3. This tension between BBN result and other 
observations might be interpreted as a 2-a evidence of the fact that further relativistic species 
are produced during the evolution of the Universe after BBN, but our viewpoint is that this 
conclusion is still statistically very weak and further data are necessary. 

|| A weaker bound is obtained by combining a larger sample of galactic HII regions, 3 He/H < (1.9 ± 0.6) • 10~ 5 



% c.L 



3 He/H > 0.95 • 10 



,-5 



(7) 
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Combining all data, under the assumption that the number of relativistic species does not 
change from BBN down to CMB and LSS formation epochs, the three neutrino scenario is still 
quite consistent, while data show a better agreement for slightly larger values Nf® = 3 — 4. 
Considering a smaller set of data would weaken the constraints on Nfp . For example, we 
mentioned that combining only the WMAP data with the recent data from the analysis of 
LRG in the SDSS [47J gives the constrain Nf = 3.7±|;f with the standard value Nf = 3.046 
well consistent with the data (see also [53J). Moreover, choosing different smaller dataset 
also affects the bound on Nf , see Table [TJ Finally, we also point out that we assumed in 
our analysis only massless neutrinos. The datasets we considered in this paper and neutrino 
oscillation experiments are indeed in agreement with a negligible (for cosmology) neutrino mass 
(see [42j and [SI]). However, considering massive neutrinos would shift our constraints to even 
higher values for Nf (see e.g. [38J). 

Apart from further data from CMB anisotropy and polarization which can be foreseen 
and available in the near future, as for example from the PLANCK experiment!^, as well as 
further studies on both Lyman-a and BAO data, we stress that new results on primordial nuclei 
abundance are needed. The large (systematic) error on the 4 He mass fraction is severely limiting 
the accuracy of any estimate of Nf at the MeV scale, while more accurate data on Deuterium, 
though always welcome, are not expected to change the overall picture in a significant manner, 
unless the accuracy in predicting the baryon density would become comparable with the one 
obtained with CMB measurements. This requires in particular, a better understanding of the 
apparent spread of individual D/H measurements in different QSO systems as due to a real non 
homogeneous distribution of primordial Deuterium in the sky or an underestimated statistical 
uncertainty in the measurements. 

Interestingly, an improved upper bound or measurement of primordial 3 He would have a 
large impact in further constraining Nf . This is unfortunately, a very difficult task, because 
of the large contamination of primordial 3 He by stellar activity. However, in the framework of 
homogeneous BBN there is a robust lower limit on this nuclide abundance, 3 He/H > 0.95- 10~ 5 , 
which is very close to the present experimental upper bound, so we can optimistically expect 
some further development in the near future. 
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